Variation of Orbital Symmetry of the Localized 3c? 1 Electron of the V 4+ Ion upon the 

Metal-Insulator Transition in V0 2 
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The electronic structure of the vanadium dioxide VO2 in the tetragonal 7? and two monocinic 
Mi and M2 structural modifications was calculated in frames of the local-density functional (LDA) 
approach and the LSDA + U formalism of correction for correlation effects. Based on the results of 
calculation, we argue in favor of the Mott-Hubbard mechanism of the metal-insulator transition, ft 
is shown that the transition is accompanied by the change of the orbital in which the 3d 1 electron 
of the V 4+ ion is localized. The symmetry of the filled orbital is determined by the local oxygen 
environment of the vanadium atoms in each structural phase. Depending on the dispersion of the 
band that corresponds to the orbital occupied, the investigated compound can be either a metal or 
an insulator. 



The metal-insulator transition in VO2 is observed at 
a temperature of 340 K and is accompanied by a struc- 
tural transition from the nonmagnetic high-temperature 
tetragonal rutile phase R to the nonmagnetic monoclinic 
Mi phase. The rutile phase contains chains of equidis- 
tant vanadium atoms along the [001] direction. The unit 
cell of the Mi phase doubles as compared to rutile, all 
vanadium atoms slightly deviate from the [001] r direc- 
tion, forming zigzags and become nonequidistant, form- 
ing pairs. 

A noncontradictory theoretical description of the 
metal-insulator transition in vanadium dioxide is still 
lacking. The origin of the nonmetallic state is explained 
now from two viewpoints. 

If one starts from the results of band calculations in the 
local-density approximation (LDA), then the Mi phase 
can be considered as a usual Peierls band insulator [ij, |2j . 
From the other viewpoint, it is impossible to ignore elec- 
tron correlations in the d shell of vanadium and they can- 
not be considered on the average, as is done in band cal- 
culations. A quarter of century ago, one more insulating 
phase (M2) with local magnetic moments at some vana- 
dium atoms was discovered. The M2 phase originates 
upon slight doping of VO2 with chromium (> 0.1 at. %) 
|3 or upon application of a small stress along the [110]# 
axis 0. In the M2 phase, in one chain, vanadium atoms 
form pairs along the [110]^ direction, as in the Mi phase 
(but without deviation from the [110],r direction), while 
in the neighboring chain the atoms of vanadium again 
become equidistant, as in the R phase, but the zigzag 
configuration is preserved, as in the Mi phase [23 . Be- 
cause of the existence of a chain of equidistant V 4+ ions, 
the band-theoretical description of the insulating state 
of the M 2 phase becomes impossible: these chains are 
magnetic (Mott-Hubbard) insulators. Therefore, theo- 
retical models were developed in which the Mi and M 2 
phases were considered as the Mott-Hubbard insulators, 
excluding the interatomic Hubbard correlations in the 
R phase At the same time, arguments and exper- 
imental evidences indicating that all structural modi- 
fications of VO2 must be considered as Mott-Hubbard 



rather than as band insulators appeared in the literature 
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Discussion of the nature of the insulating state in VO2 
continues with allowance for new experimental and the- 
oretical data on the electronic structure of the R and 
Mi phases. From comparison of the results of spectro- 
scopic measurements and band calculations, it was con- 
cluded in that the electronic structure of VO2 is well 
described within the local-density approximation with- 
out an account of any correlation effects. Experimental 
measurements of the temperature dependence of the ve- 
locity of surface acoustic waves in the VO2 single crys- 
tal near the metal-insulator transition permitted Mau- 
rer et al. to conclude that it is unlikely that the 

driving force of the metal-insulator transition arc elec- 
tron correlations. Measurements of the delay time in the 
VO2 metal/oxide/metal structures and the investigation 
of the influence of the injection of carriers on switching 
and metal-insulator transition in planar switching devices 
fl4| also confirm exclusively the Mott-Hubbard mecha- 
nism of transition in VO2. It is known that the usual 
band calculations indicate only a tendency for the forma- 
tion of an energy gap in the Mi phase of VO2 Q and give 
a resulting metallic spectrum in contrast to the experi- 
ment. An attempt to resolve this contradiction within 
the band approach was made. Thus, the introduction of 
the electron-phonon interaction into the combination of 
the three-dimensional periodic-shell model and the Xa 
cluster method of discrete variation resulted in the ap- 
pearance of a gap in the band spectrum [l5| . The intro- 
duction of the correction for self-energy within the model 
GW scheme also satisfactorily reproduces the properties 
of the non metallic Mi phase of VO2 observed experi- 
mentally 0]. Note that no theoretical calculations of 
the electronic structure of the insulating M2 phase are 
available at present. 

In this work, we tried to analyze in detail the electronic 
structure of the R, Mi and M 2 phases of VO2 (for the 
first time in the case of the M2 phase), in each case paying 
attention to the orbital that is filled by the sole electron 
of the V 4+ ion. With this aim in mind, we used the first- 
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FIG. 1: The crystal structure of the R phase of VO2. Large 
spheres are the vanadium atoms, small spheres are the oxygen 
atoms. The axes of the crystallographic coordinate system are 
labeled by capital letters. The axes of the local systems of 
coordinates for the corner and body-center vanadium atoms 
are labeled by small letters. 



principles linearized MT-orbital method of calculation in 
the tight-binding approximation (TB LMTO) [H [ufy . 
The nonmagnetic calculations were carried out within 
the conventional formalism of the local-density functional 
(LDA). In magnetic calculations, the electron correla- 
tions were taken into account within the LSDA+U ap- 
proach 0, |2(J ■ We used the experimentally determined 
values of the lattice parameters available in the literature 
(from 0| for the R phase, from |22| for the Mi phase, 
and from for the M2 phase). In calculations, the radii 
of the atomic spheres of V and O were chosen to be equal 
with a 16% overlap; empty spheres were introduced to fill 
the volumes of the atomic spheres to the volume of the 
unit cell; the basic functions were 4s, Ap, and 3d for vana- 
dium, 3s, 2p, and 3d for oxygen, and Is, 2p, and 3d for 
the empty spheres. Within the LSDA formalism, we cal- 
culated the Coulomb U and exchange J interactions in 
the d shell of vanadium. The value of U was calculated 
with allowance for screening of the ti g electrons by the 
e g electrons j23|. The values U — 3.8 eV and J = 1.0 eV 
obtained were used unchanged in calculations of the elec- 
tronic structure of all phases j^. Thus, our calculations 
do not contain adjustable parameters and, in this sense, 
are the first-principles. 

The unit cell of the high-temperature R phase of VO2 
(space group PA2 / mnm) contains two formula units. The 
vanadium atoms form a body-centered tetragonal lattice 
and are surrounded by slightly deformed oxygen octahe- 
dra. In such octahedra, the V-0 distances in the direc- 
tion normal to the tetragonal c axis are slightly greater 
than the V-0 distances in the plane parallel to the c axis 
(Fig. 2] and table). In the same plane, the O-V-0 angles 
are 84.2° and 95.8°, while the other O-V-0 angles are 



equal to 90°. The octahedra surrounding the vanadium 
atoms in the vertices and in the center of the unit cell are 
rotated by 90° about the c axis. In what follows, we will 
present the results for each vanadium atom in the local 
coordinate system, in which the z axis is directed toward 
one of the most remote oxygen atoms (apical atom) in 
the octahedra. The direction of the chain of vanadium 
atoms along the c axis is the bisector of the angle which 
is made by the x and y axes of the local coordinate sys- 
tem (see Fig.nj. The Euler angles used for passage from 
the crystallographic coordinate system to the local one 
for the vanadium atom at the position (0; 0; 0) are the 
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Let us clarify the terms that will be used in the inter- 
pretation of the results. By occupancy of the ith orbital, 
we mean the integral iVj = J^ F rii(e)de where rii(e) is 
the corresponding partial density of states, and Ep is the 
Fermi energy. By the center of gravity of the ith band, we 
mean the ratio e^ 9 =< Ej >= e-ni{e)de/ m(e)de, 
where E\ and £2 are the boundaries of the energy interval 
under consideration. By the dispersion of the ith band, 
the expression u, = y/ < ef > — < £j > 2 is meant. Be- 
cause of the distortion of the oxygen octahedra, the ti g 
orbitals of the vanadium atom become nondegenerate. 
The xy, yz — zx, and yz + zx wave functions, where x, 
y, and z represent the basis of the local system of coor- 
dinates, will be the natural basis for the tig orbitals. 

The results of the nonmagnetic LDA calculation are 
shown in Fig.|3 Our LDA calculations show: (1) owing 
to the elongation of the oxygen octahedra along the local 
z axis, the yz — zx orbital is most populated; (2) the 
center of gravity of the yz — zx band lies most deeply 
in the energy range from -8 to +6 eV; (3) in the energy 
interval from -1.5 to +1.8 eV, the dispersion of the yz—zx 
band is approximately one and half times greater than 
that for the xy and yz + zx bands. 

We also carried out LSDA+U calculations for the R 
phase of VO2 for the ferromagnetic and antiferromag- 
netic types of ordering of magnetic moments of vanadium 
ions in the chains located along the tetragonal c axis. In 
the first case, the half-metallic ferromagnetic state is ob- 



TABLE I: The V-O distances (A) for the R and Mi phases 
and the average distances for the M2 phase of VO2. The 
"apical" and "in-plane" oxygen atoms are considered in the 
local coordinate systems: the z axis passes through the apical 
oxygen and the x and y axes pass through the in-plane oxygen 
atoms. 
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FIG. 2: The total (per formula unit) and partial (per atom) 
densities of states for the R phase of VO2 (left column) and 
the ti a orbital partial density of states for the vanadium atoms 
(right column) . The positions of the Fermi level are shown by 
dotted lines. Arrows in the right column indicate the centers 
of gravity of the orbital partial densities of states in the energy 
range from -8 to +6 eV. 

tained with the spin magnetic moment of the d shell of 
vanadium equal to /i=0.99 /xg. In the second case the 
antiferromagnetic state was found to be insulating with 
an energy gap of 0.38 eV and ^=0.90 [Ib- The second 
solution was lower than the first one in the total energy. 
For both solutions the sole d electron of the V 4+ ion was 
localized in the yz — zx orbital (Fig. [3J|. This naturally 
follows from the above results ((1) and (2)) of the LDA 
calculations: within the LSDA+U approach, the orbital 
that already was most populated and whose center of 
gravity was already the lowest among all the ti g orbitals 
is occupied. The third result (the maximum dispersion of 
the yz — zx band) has preserved the ferromagnetic state 
as the metallic one even within the LSDA+U calcula- 
tion and led to the appearance of the antiferromagnetic 
insulating state with a small energy gap|30j. 

The fact that the LSDA+U calculations predict the an- 
tiferromagnetic insulating state for the R phase of VO2 
instead of the metallic one reflects the known tendency 
of the LSDA+U method to overestimate the electron lo- 
calization as a result of consideration of the on-site spin- 
and orbit-dependent Coulomb repulsion in terms of the 
mean-field theory. At the same time, the LSDA+U ap- 
proach gives an accurate description of the orbital nature 
of the localized electrons |3l| . Our calculations of the elec- 
tron structure of VO2 in the R phase predict localization 
of the d electron of the V 4+ ion on the yz — zx orbital. 
As will be seen below, this result is a key one for the 
description of the metal-insulator transition in VO2. 

Analysis of the partial densities of the ti g states of 
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FIG. 3: Ordering of the occupied yz — zx orbitals in the 
chain of the vanadium atoms obtained as a result of the 
LSDA+U calculation for the R phase of VO2 For designa- 
tions, see Fig. The orbitals shown represent the angular 
distribution of the d electron density of the vanadium atoms 
p{S,4>) = J2 m , m ' Qm,m'Y m {9,(f))Y m ,(9,(l>), where Q m>m > = 
n T , — n* , is the difference between the population matri- 
ces of the d states of the vanadium atoms with spins up and 
down obtained by self-consistent calculation, and Y m (9, </>) are 
the corresponding spherical harmonics. 



vanadium, obtained in the nonmagnetic LDA calcula- 
tion for the Mi phase and shown in Fig. indicates 
a different sequence of filling orbitals, as compared to 
the R phase. The M\ phase (the space group P2\jc) is 
characterized by a pairing of the vanadium atoms along 
the [0.842; 0; —0.539] vector in the crystallographic co- 
ordinate system: there are short and long vanadium- 
vanadium bonds. Moreover, each pair of vanadium atoms 
is inclined relative to this vector in the plane with the nor- 
mal [0.381; 0.708; 0.595] for one chain of vanadium atoms 
and to the plane with the normal [0.381; -0.708; 0.595] 
for the neighboring chain. The difference in the oxygen 
surroundings of the vanadium atoms is not only in the 
variation of the O-V-0 angles, but also in that the V- 
O distance to the apical oxygen atoms in the M\ phase 
became smaller than to the in-plane oxygens (see table) . 
We define the local system of coordinates for the M\ 
phase in the following way: the z axis is perpendicular 
to the plane of the vanadium atoms that belong to the 
same chain, and the x and y axes are deflected by an 
angle of 45° from the direction of the chain (Fig. |SJ. In 
other words, the local system of coordinates for the Mi 
phase is chosen similar to that for the R phase. Then, 
for the vanadium atom in the (0.259; —0.025; —0.136) po- 
sition, the Euler angles for passing from the crystallo- 
graphic system of coordinates to the local one will be 
(0.6587r;-0.2977r;-0.9847r). 

Comparing the partial densities of the ti g states of 
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FIG. 4: The total (per formula unit) and partial (per atom) 
densities of states for the M\ phase of VO2 (left column) 
and the t2 g partial densities of states for the vanadium atoms 
(right column). The Fermi level is shown by dotted lines. 

vanadium for the R and Mi phases (see Figs. and |3J, 
note a substantial difference in the density of states for 
the xy orbital. Owing to dimenzation of the vanadium 
atoms, the corresponding band for the Mi phase is split 
into bonding and antibonding subbands and, which is 
more important, the filling of the xy orbital becomes 
maximum between the t2 g orbitals. The nonmagnetic 
LDA calculation does not reproduce the insulating nature 
of the energy spectrum for the Mi phase of VO2. Yet, it 
is clear that the allowance for correlations in calculations 
will not change this key result, i.e., the preferential filling 
of the orbitals of the xy symmetry. 

The conclusion concerning filling of the 3d 1 xy states 
of the Mi phase is supported by two factors. First, the 
distortions of the oxygen octahedra around the vana- 
dium atom occur in such a way that the maximum V- 
O distances are observed in the xy plane in the local 
system of coordinates (see table), which, in turn, low- 
ers the energy of the xy orbital. Second, we performed 
the LSDA+U calculations for the Mi structure to de- 
termine which of the orbitals will be filled by the local- 
ized electron |3^. Note that the LSDA+U calculation 
predicts for the Mi phase the antiferromagnetic ground 
state with a gap equal to 0.55 eV and the spin magnetic 
moment of the d shell of vanadium /z=0.81 \xb instead of 
the nonmagnetic insulating state with a gap of 0.2-0.7 eV 
observed experimentally [9, 27]. Quantum corrections, 
probably, could suppress the long-range magnetic order 
obtained in the LSDA+U calculation, and the calculated 
value of the gap is within the limits of the experimental 
one. Yet, turning to the discussion of the orbital nature 
of the localized d electron, one can see from Fig. 5 that 
the xy orbital became localized. Both factors permit us 




FIG. 5: Ordering of the occupied xy orbitals in the chain of 
vanadium atoms obtained as a result of the LDA+U calcula- 
tion for the Mi phase of VO2 . The axes of the cristallographyc 
system of coordinates are marked by capital letters. The axes 
of the local system of coordinates for the vanadium atoms are 
marked by small letters. 

to conclude that upon transition of VO2 from the R to 
the Mi phase, along with a change in the ratio of the 
distances from the vanadium atom to the in-plane oxy- 
gen atom and to the apical oxygen atom, an important 
transformation of orbital symmetry of the localized 3d 1 
electron of the V 4+ ion occurs, namely, from the yz — zx 
to the xy. 

We believe that the Peierls mechanism is not respon- 
sible for the formation of the insulating state of the Mi 
phase. Even the LDA results show that the most filled 
orbitals in both cases are different: the xy orbital for the 
Mi phase and the yz — zx orbital for the R phase. It is 
known (for example, from the theoretical investigations 
of Cr02 [3 which also has a crystal structure of the ru- 
tile type) that the yz — zx band has substantial dispersion 
and is strongly hybridized with the 02p states, forming 
the 7T bond with them. Therefore, it is weakly localized 
even in the presence of strong correlations and the energy 
spectrum may remain metallic in the case of the prefer- 
ential filling of this band. At the same time, the xy band 
is narrower than the yz — zx one and more sensitive to 
correlations. This results in localization of electrons on 
the xy orbital, its complete filling, and the appearance 
of a gap in the energy spectrum. It is what we obtained 
in our LSDA+U calculations. Thus, we can conclude 
that it is the Mott-Hubbard mechanism that is respon- 
sible for the formation of the energy gap in VO2. The 
Peierls's mechanism can be used only for explanation of 
the insulating state of the Mi phase, but it is completely 
unacceptable for the M2 phase, in which both dimerized 
and undimerized chains of vanadium atoms exist simul- 
taneously. Taking into account which orbitals arc filled 
by the d electron in each particular case and which are 
the properties of these orbitals, the Mott-Hubbard mech- 
anism can be applied for the description of both metallic 



and insulating states of all crystal phases of VO2. 

The calculations of the electronic structure of the Mi 
phase (space group C2/m) fully confirm our assumption 
about the mechanism of formation of energy bands in 
VO2. Both dimerized chains of vanadium atoms (similar 
to the chains in the M\ phase) and the undimerized ones 
(as in the R phase) exist in the crystal structure of the 
Mi phase. As in the Mi phase, but contrary to the 
R phase, the V-0 distances to the apical oxygen atoms 
in the oxygen octahedra which surround the vanadium 
atoms in both types of chains are smaller than to the in- 
plane oxygens. The Euler angles for the passage to the 
local system of coordinates for the vanadium atoms in 
the dimerized and undimerized chains (Fig. 6) are given 
by the sets of (-ir; -f ; f ) and (0; -f ; f ). 

The "nonmagnetic" calculations in the local-density 
approximation give the metallic ground state, which con- 
tradicts the experiment. Owing to a substantial distor- 
tion of the crystal structure and, correspondingly, com- 
plicated hybridization of the d states of vanadium with 
each other, it is impossible here to consistently separate 
the partial densities of the ti g states in the cubic harmon- 
ics, which we had done, for example, for the R phase. The 
LSDA+U calculation is the sole possibility of determin- 
ing which localized ti g orbital will be filled. 

The results of such a calculation are shown in Fig. 
They demonstrate the occupied orbitals of vanadium that 
in the Mi phase are the xy orbitals both in the dimer- 
ized and undimerized chains. According to calculations, 
in the case of antiferromagnetic ordering the spin mag- 
netic moment of the d shell of vanadium that belongs 
to the undimerized chain is equal to fj, = 0.89 fj,B and 
the arising energy gap is 0.93 eV. We also calculated 
the exchange-correlation parameter J of the Heisenberg 
model for the vanadium atoms in the undimerized chain 
(the scheme of calculation within the LSDA+U method 
is described in j^jj. Its value was found to be 205 K and 
the sign corresponds to the antiferromagnetic exchange. 
The value obtained is close to the result of Q where, 
from an analysis of nuclear magnetic resonance experi- 
ments, it is shown that the V 4+ ions in the undimerized 
chains behave as Heisenberg chains with spin S = 1/2 
and antiferromagnetic exchange J«300 K. 

Thus, the principal difference between the rutile R 
phase and the monoclinic Mi and M2 phases is which of 
the localized ti g orbitals is filled by the sole 3d electron 
of the V 4+ ion: yz — zx or xy. The choice of the filled 
orbital is determined by the local oxygen surroundings 
of the vanadium atoms. Since these bands have different 
properties relative to localization (or dispersion), the cor- 
relation effects variously influence the energy spectrum 
(metal or insulator) depending on which band is filled. 
It is correlations (the Mott-Hubbard mechanism) that 
are responsible for the insulating properties of the M\ 
and Mi phases of VO2. Correlations are also present in 
the metallic R phase. Owing to the localization of the d 




FIG. 6: Ordering of the occupied xy orbitals in the dimerized 
(upper left part of the figure) and the undimerized (lower right 
part of the figure) chains of the vanadium atoms obtained by 
the LSDA+U calculation for the M2 phase of VO2. The axes 
of the crystallographic system of coordinates are marked by 
capital letters. The axes of the local systems of coordinates 
are shown by small letters. 



electron in the yz — zx orbital, they do not lead to the 
appearance of a gap in the band spectrum. 
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The diagram of mutual arrangement of V atoms in the 
three VO2 phases can be found, for example, in [[f|, Fig. 
215]. The phase diagram is shown in [0, Fig. 217]. The 
triclinic T phase of VO2 will not be considered in this 
work. 

The existing experimental estimates of the Coulomb in- 
teraction U in VO2 give 2.1 eV for the insulating 
phase and 1.3 eV for the metallic phase. These values 
are slightly below that used in calculations U = 3.8 eV. 
For U < 3.0 eV, the antiferromagnetic state of the R 
phase in the LSDA+U calculation is metallic. See also 
[29]. 

This is shown in the series of our works devoted to the 
investigation of orbital filling in Cr02 24] KCUF3 per- 
ovskites [20|. PrMnOs |25l| . layered vanadates of the 
CaV„02n+i-type [2(|, etc. 

In calculating the Mi and M2 phases of VO2, it is incor- 
rect to substantiate conclusions by comparing the centers 
of gravity of orbitals and the dispersions of the corre- 
sponding bands obtained in nonmagnetic LDA calcula- 
tions in view of the lower symmetry of these phases and 
substantial hybridization of the states under investiga- 
tion. 



